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1.  INTRODUCTION 


In  the  mid  1950s  four  shell  types,  exhibited  unusual  behavior  that  involved  the  movement  of  their 
interior  parts  (Boyer  1955a,  1955b;  Roecker  and  Boyer  1957;  Karpov  and  Bradley  1958).  In  1977  a 
simple  theory  was  developed  to  relate  the  motion  of  internal  parts  to  the  angular  motion  of  the  spinning 
projectile  (Murphy  1977).  This  theory  assumed  a  circular  lateral  motion  of  the  internal  component 
combined  with  a  coning  motion  of  its  axis  of  symmetry.  Both  of  these  motions  are  at  a  coning  frequency 
of  the  spinning  shell.  The  frequency  of  the  internal  motion  was  assumed  to  be  the  fast  frequency  of  the 
projectile’s  free  angular  motion,  the  amplitude  was  set  by  the  availaUe  tolerances,  and  the  phase  of  the 
motion  was  induced  by  internal  friction. 

The  theory  predicted  a  destabilizing  side  moment  and  a  spin-up  moment  equal  to  the  product  of  the 
side  moment  and  the  magnitude  of  the  projectile’s  angular  motion.  Later  it  was  shown  that  this  relation 
between  side  moment  and  spin  moment  is  valid  for  any  steady  motion  of  an  internal  payload  (Murphy 
1989).  For  60°  phase  angles,  all  observed  instabilities  could  be  explained  by  the  measured  clearances  and 
the  measured  despins  of  the  unstable  projectiles  were  in  good  quantitative  agreement  with  the  theory’s 
prediction. 

The  theory,  however,  made  no  attempt  to  predict  the  amplitude  and  phase  of  the  internal  motion.  For 
a  mass  suf^rted  by  an  elastic  internal  beam,  W.  R.  Chadwick  (1975)  did  develop  a  complete  theory  for 
the  combined  motion  of  the  internal  mass  and  the  projectile.  Unfortunately,  this  work  was  marred  by  two 
errors  in  the  proper  treatment  of  beam  damping,  and  his  equations  did  not  predict  flight  results. 

In  this  report,*  we  develop  the  correct  equations  of  motion  for  the  projectile  and  its  internal,  beam- 
supported  mass.  Using  Chadwick’s  influence  coefficients,  we  compute  the  elastic  response  of  various 
beams  with  a  proper  expression  for  beam  damping.  The  beams  considered  are  forward-facing  and 
rearward-facing  cantilever  beams,  fixed-pirmed  beams,  and  fixed-fixed  beams. 

The  theory  and  associated  computer  runs  show  that  the  instabilities  induced  by  beam-moimted  masses 
are  strongly  controlled  by  beam  elastic  characteristics  and,  very  importantly,  the  form  and  amplitude  of 


*  A  shorter  version  of  this  report  has  appeared  as  an  AIAA  preprint  (Murphy  1992). 
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the  beam  damping.  Small  beam  damping  has  little  effect  on  projectile  flight  stability,  but  moderate 
damping  can  cause  instability  when  the  beam  is  sufficiently  soft 

The  inertia  properties  of  the  internal  mass  determine  the  relative  importance  of  beam  deflection  or 
beam  cant  on  projectile  stability.  Theory  indicates  design  tradeoffs  for  beam  characteristics  and  mass 
inertial  properties. 

2.  EQUATIONS  OF  MOTION 


In  Murphy  (1977),  the  internal  component  is  assumed  to  perform  a  known  lateral  and  angular  motion. 
In  complex  notation,  the  lateral  motion  is  described  in  nonspinning  coordinates  by 

Yc  +  K  =  (1) 

and  the  angular  motion  is  described  by  the  component’s  axis  of  symmetry  orientation  angles  with  respect 
to  the  projectile’s  axis  of  symmetry: 

Ty  +  iYz  =  r.  (2) 

The  projectile’s  motion  relative  to  its  flight  path  is  given  by  its  angles  of  attack  and  sideslip  in 
nonroUing  coordinates: 

p  +  ia  =  |.  (3) 

These  geometrical  quantities  are  shown  in  Figures  1  and  2.  The  usual  linear  analysis  yields  the  following 
equation  for  the  projectile’s  motion: 


Itf  -  (Aq  +  iL,)|  -  (A„  +  iAp„)|  =  B,  E  I,,  r  -  ip,  I,,  f. 


(4) 


The  angular  motion  of  a  solid  projectile  (E  =  F  =  0)  is  usually  described  as  the  sum  of  two  coning 
motions  of  the  form 


I  =  I 

j  =  l,2 


KjC 


(5) 


21, 


1  ± 


f 


l/s„ 


(6) 


2 


3 


(7) 


and 


L 


2 

X 


41, A„ 


_L  =  X;  = 


(8) 


Since  the  fast  mode  is  usually  the  one  adversely  affected  by  payload  motion,  a  simple  one-mode  motion 
is  assumed,  and  the  internal  component  is  assumed  to  perform  constant-amplitude  motion  at  the  fast 
frequency. 


and 


(9) 


E  = 


(10) 


When  equations  9-10  are  substituted  in  equation  4,  it  can  be  shown  that  the  firequency  and  damping 
of  the  fast  mode  are  changed  by  the  following  amounts; 


and 


where 


Ki(2I,(^l  -LJ 


AXj 


K,(2I,<i.i  -Lj’ 


(11) 


(12) 


^1  “  ycos  <j).y  -  Bg  e  cos  4>e 
81  =  8^7  sin  -  Bg  (j»  j  e  sin  0g 

~  ^xc  Pc  ~  ltc9  1 
Be  = 
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y  and  e  and  their  phase  angles  were  estimated  firom  available  clearances  and  good  engineering  predictions 
were  obtained.  In  this  report,  we  will  predict  these  quantities  for  a  component  mounted  on  an  elastic 
beam. 

3.  COMPONENT  EQUATIONS  OF  MOTION 

The  lateral  motion  of  an  internal  component  can  be  described  in  a  nonrolling  coordinate  system  whose 
X-axis  is  aligned  along  the  projectile’s  axis  of  symmetry.  The  angular  velocity  of  these  coordinates  can 
be  specified  by 

Q  =  0  +  i\jr. 


The  lateral  motion  of  the  component  can  be  easily  expressed  in  terms  of  the  transverse  components  of  the 
force  and  moment  exerted  by  the  beam  on  ihe  component:* 


and 


“ 

r 

“ 

nic 

{E  +  V 

1  -  iQ 

“  iXcQ 

_ 

1 

- 

(13) 


Itc(l'~*Q)“^Pclxc(^  iQ)=  (^C  ^^zc)' 


(14) 


The  elastic  part  of  the  beam  force  and  moment  can  be  computed  by  simple  beam  theory  to  be  a  linear 
function  of  E  and  T.  For  the  coning  motion  of  equations  9-10, 


arid 


[Pyc  +  ^^zc]  elastic 


b,,  {Beil’s  +  bi2ye‘''*^ 


[Myc  +  elastic 


i  [b2i{eei^s  +  b22Yei^Y]e*‘’’‘, 


(15) 


(16) 


*  See  Appendix  A  for  derivation  of  Equations  13-14.  Appendix  B  outlines  the  derivation  of  Table  1. 
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where  bi2  -  and  the  linear  coefficients,  bjj ,  can  be  computed  for  a  variety  of  beam  supports  by 
taking  the  matrix  inverse  of  the  influence  coefficients,  ajj,  given  in  Table  1. 

The  effect  of  beam  damping  can  be  approximated  by  linear  terms  in  the  derivatives  of  E  and  F.  Since 
the  beam  is  spinning  with  the  projectile,  these  derivatives  must  be  taken  in  projectile-fixed  coordinates 

Table  1.  Influence  Coefficients  for  Various  Beams 


Fixed-fixed 


Fixed-pinned 


Pinned-fixed 


Fixed-free 


Free-fixed 


Pinned-pinned 


3(a  +  b) 


a^b^(3a  +  4b) 


a2b^(4a  +  3b) 


a2b(2b2  -  a^) 


ab^Cb^  -  2a ^) 


a(a^  +  4b  ^) 


b(4a^  +  b^) 


3(a  +  b) 


ab(b  -  a) 


3(a  +  b) 


a^  +  b^  -  ab 


3(a  +  b) 


that  are  spinmng  with  respect  to  earth-fixed  axes.  A  coning  frequency  j  in  the  nonspinning  system 
would  appear  in  a  spinning  system  as  -  p  where  p  is  the  projectile  spire  The  beam  damping  can  be 
inserted  in  equation  (15)  by  replacing  by  by  bjj  j^l  +  i  dg  2  -  p)  p  ”  and  in  equation  (16)  by  replacing 
*>22  [l  +  idy  (^i  -  p)p"^j. 

If  the  small  effect  of  the  lift  force  is  neglected. 


I  =  iQ. 
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Equations  13-16  can  be  combined  to  yield 


+  bij[l  +  idgC^i  -  P)P“^]  +  bi2Ye‘'>’Y  =  mj.Xc((^i)2Ki  (18) 


and 


b2i«ee^*>’e  +  -  01  By  +  b22  [l  +  idy  (4>i  -  p)  P~^] 
=  -  i  01  By  Ki. 


(19) 


When  beam  parameters,  component  inertia  properties,  and  coning  motion  parameters  are  inserted  in 
equations  18-19,  amplitude  and  phase  angles  of  the  component  motion  can  be  determined.  From  these, 
the  stability  of  the  coning  motion  follows  from  equation  12. 

4.  DISCUSSION 

To  illustrate  these  results,  we  will  consider  an  8-in  projectile  with  a  particular  rotationaUy  symmetric 
internal  component.  The  appropriate  physical  and  aerodynamic  properties  of  this  projectile  and  its 

component  are  given  in  Table  2.  The  nominal  value  of  the  beam’s  stiffness  is  given  as  El*  = 
59,800  Ib/ft^.  In  Figures  3  and  4,  the  projectile  damping  is  plotted  as  a  function  of  beam  elasticity  for 
all  six  possible  beam  types.  When  the  beam  stiffness  is  reduced  to  10%  of  the  Table  2  value,  all  beam 
types  can  cause  flight  instability. 


Table  2.  Test  Case  Values 


In  Figure  5,  the  size  of  the  component  motion  is  plotted  as  a  function  of  beam  stiffiiess  for  the  fixed- 
fixed  beam.  For  small  coning  motion  (Kj  <  0.1),  this  motion  is  quite  small  and  only  die  relatively  large 
value  assigned  to  causes  this  motion  to  induce  flight  instability. 


Figure  5.  Comtxjnent  cant  and  deflection  as  functions  of  beam  stiffness  for  a  fixed-fixed  beam. 

Figure  6  shows  die  beam  damping  effect  on  flight  stability  for  the  fixed-fixed  beam  with  EI/EI  *  = 
0.1.  The  nominal  value  of  d  =  0.5  is  shown  for  the  equal  damping  curve  d^  =  dy  =  d.  The  odier  two 
curves  consider  the  solo  effects  of  continued  damping  (dy)  and  of  deflection  damping  (d^),  respectively. 
Qearly,  cant  damping  can  have  a  greater  adverse  effect,  and  0.5  is  the  worst  value  for  dg  =  dy. 


Figure  6.  Projectile  damning  as  a  function  of  beam  damping  for  a  fixed-fixed  beam  and  three  types  of 
damping:  EI/EI*  =  0.1. 
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Intuitively,  one  might  think  the  instabilities  shown  by  Figures  3  and  4  are  due  to  reducing  the  natural 
ftequency  of  the  beam  to  the  coning  frequency  of  the  projectile.  In  Figure  7,  the  natural  beam  frequency, 
of  the  three  beam  types  of  Figure  3  is  computed,  and  the  projectile  damping  is  plotted  as  a  function 

of  (j)  For  the  fixed-fixed  beam,  instability  occurs  when  the  coning  frequency  is  only  1 0%  of  the  beam 
frequency! 


Figure?.  Pro]ectile  damping  as  function  of  coning  freauencv  for  three  nonninned  beam  tvnes. 
5.  CONCLUSIONS 

(1)  The  steady-state  motion  of  internal  mass  on  an  elastic  beam  can  be  computed. 

(2)  For  appropriate  values  of  beam  damping,  flight  instabilities  can  be  predicted. 

(3)  These  instabilities  can  occur  for  frequencies  less  than  10%  of  the  beam  natural  frequency. 
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APPENDIX  A: 
EQUATIONS  OF  MOTION 


13 


Intentionally  left  blank. 


14 


The  complete  projectile  is  assumed  to  consist  of  a  body  with  mass  m^  and  an  internal  component  of 
mass  m^.  mounted  on  an  elastic  massless  beam.  The  total  mass  of  the  projectile  is  m  =  m^.  +  m^.  The 
center  of  mass  (CM)  of  the  internal  component  is  located  a  distance  x^.  from  the  CM  of  the  projectile,  and 
the  CM  of  the  body  is  a  distance  x,,  from  the  CM  of  the  projectile  (m^  x^  +  m^.  x^  =  0).  If  is  the 
transverse  moment  of  inertia  of  the  component  relative  to  its  CM  and  Ijj,  is  the  transverse  moment  of 
inertia  of  the  body  relative  to  its  CM,  then  Ij,  the  transverse  moment  of  inertia  of  the  projectile  relative 
to  its  CM,  is  given  by 


It  =  Itb 


+  m,,x^ 


+  I 


tc 


mcX^ 


(A-1) 


The  projectile’s  axis  of  symmetry  can  be  located  relative  to  fixed  axes  by  its  pitch  and  yaw  angles, 
0  and  'F,  and  its  orientation  relative  to  the  velocity  vector  can  be  defined  by  its  angles  of  attack  and  side¬ 
slip,*  a  and  p.  The  internal  component  is  assumed  to  be  a  cylinder  or  disk  whose  axis  is  located 
relative  to  the  projectile’s  axis  by  the  angles  jy  and  The  transverse  displacement  of  the  component  CM 
relative  to  the  body’s  axis  of  symmetry  is  specified  by  the  displacements  and  z^.  All  these  geometrical 
quantities  are  shown  in  Rgures  1  and  2. 

The  projectile  has  lateral  aerodynamic  forces  and  moments  acting  on  it,  represented  by  Fy,  and  My, 
Mj,  respectively.  The  elastic  beam  exerts  forces  and  moments  on  the  internal  component  when  it  is 
deflected;  these  are  represented  by  Fy^,  F^^,  and  My^,  M^,.,  respectively.  For  simplicity,  aerodynamic  drag 
is  neglected  and  all  angles  are  assumed  to  be  smaU. 

The  equations  of  motion  of  the  projectile  body  and  internal  component  are  as  follows: 


mb[v  (p  +  rj/) 


(A-2) 


*  The  tilde  superscripts  are  used  to  indicate  that  these  angles  are  measured  in  nonspinning  coordinates  and  not  the  usual  missOe- 
fixed  coordinates. 


15 


nib  [V  (a  -  e)  -  Xb0]  =  (A-3) 

Itb  6  +  Pb  Ixb  ¥  =  My  -  Myc  +  Xb  F^  +  (x^  -  Xb)  F^^ ,  (A-4) 

^tb  ¥  -  Pb  Ixb  6  =  M^  -  -  Xb  Fy  -  (x^  -  Xb)  Fy^ ,  (A-5) 

“c  [y’c  +  V (p  +  v)  +  x^  v]  =  Fyc ,  (A-6) 

nic  [z’c  +  V(d  -  e)  -  x^  0]  =  F^^,  (A-7) 

Itc  (6  -  y'z)  +  Pc  Ixc  (¥  +  Yy)  =  My,,  (A-8) 

and 

Itc  (¥  +  Yy)  -  Pc  Ixc  (0  -  Yz)  =  M,,  ,  (A-9) 


where  I^b  ^xc  die  axial  moments  of  inertia  of  the  body  and  component,  respectively. 

The  transverse  force  on  the  component  can  be  eliminated  between  equations  (A-2)-{A-3)  and 
(A-6)-(A-7),  and  complex  variables  can  be  introduced  to  yield  the  following  force  equation  for  the 
projectile: 

mV  (|  -  i  q)  =  Fy  +  i  F^  -  m,«E,  (A-10) 
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where 


and 


I  =  P  +  ia, 

Q  =  9  +  ixjr, 

E  =  (yc  +  iZc) 


The  external  aerodynamic  force  can  be  approximated  by  the  linear  normal  force,  and  the  small  component 
inertia  force  can  be  neglected  to  yield  the  simple  relation 


mV 


(i-iQ). 


(A-ll) 


where 

Fn  =  (i/2)pV2sCn„. 

Equations  (A-2)-(A-9)  can  now  be  combined  to  eliminate  the  transverse  forces  and  moments  acting 
on  the  internal  component  The  resulting  moment  equation  for  the  projectile  is 


Ij  Q  -  iL^Q  =  My  +  iM^ 


iBgE  -  ilt^.r  - 


Pcixc 


r. 


(A-12) 


where 

“  Pb  ^c’ 


Be  =  Xc  C, 
and 

r  =  Yy  +  iYz- 
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The  external  aerodynamic  moment  can  be  replaced  by  the  usual  linear  moment  expansion  and  Q  can  be 
eliminated  by  the  use  of  equation  (A- 11)  to  yield  the  following  form  of  the  projectile  momait  equation; 


It  f  -  K  -  iL,)  i  -  (a„  .  i  ApJ  I  -  B,  E  -  I,,  f  +  i  p^  r 


0, 


(A-13) 


where 


A-  =  (l/2)pS«2  V 


^Ma 


a 


A„  =  (l/2)pS«v2C^,^, 


and 


Ap„  =  (1/2)  p  S  V  p, 


Cm  +  k,  C»j 
Mp^  a  iNct 


The  complex  equations  for  the  motion  of  the  internal  component  can  be  written  from 
equations  (A-6)-(A-9)  as  follows: 


f 

me 

«  E  +  V 

i-iQ 

-  iXeQ 

_ 

and 


(A-14) 


Itc  (f  -  iQ)  -  iPcIxc  (r  -  iQ)  =  -  i  (My^  +  i  M^^).  (A-15) 


The  motion  of  the  internal  component  is  determined  by  equations  (A- 11),  (A-14),  and  (A-15)  when  the 
component  forces  and  moments  are  specified  and  the  projectile  motion  is  known. 
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APPENDIX  B: 

BEAM  INFLUENCE  COEFHCIENTS 
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Intentionally  left  blank. 
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If  a  mass  is  moxmted  on  an  elastic  beam  and  the  beam  is  deformed,  the  beam  exerts  a  force  and 
moment  on  the  mass.  The  components  of  this  force  and  moment  a^jear  in  equations  (A-14)  and  (A-15) 
and  it  is  necessary  to  express  these  quantities  as  functions  of  the  displacement  (£y,  e^,)  and  rotation  (Yy, 
Yj)  of  the  mass.  For  an  elastic  beam,  these  are  linear  relations: 


F,c*‘F«-l>u<E*b,2r 


(B-1) 


and 


Myc  +  »Mzc  =  i(b2ifiE  +  b22r). 


(B-2) 


where  bj2  =  ^21. 


For  simple  beam  theory,  it  is  easy  to  express  the  displacement  and  rotation  as  functions  of  applied 
force  and  moment  (F^,  M^).  In  the  xy  plane,  this  yields  the  following  linear  relations  in  terms  of  the  beam 
influence  coefficients: 


=  ^11  Fa  +  ai2  (B-3) 

and 

Yy  =  a2lFa  +  322  M^. 

Fj  and  are  the  negatives  of  the  corresponding  force  and  moment  on  the  internal  component  (Fy^, 

Thus,  by  equations  (B-1)  to  (B-4),  the  b^’s  are  related  to  tiie  ajj’s  as  follows: 


toll 


^22  - 


and 


D  -  a22  3i2  “  (^12)^- 


(B-5) 

(B-6) 

(B-7) 

(B-8) 


21 


We  will  assume  the  beam  has  a  length  a  +  b;  the  internal  component  and  its  associated  applied  force  and 
moment  are  located  a  distance  a  from  the  forward  end  of  the  beam  and  a  distance  b  from  its  rear  end. 

The  influence  coefficients  can  be  computed  for  various  beams  by  means  of  the  weU-known  elastic 
beam  equation: 


where 


El 


(B-9) 


and 


^  ^y  0  ^  0 

o<x<  a 

-xFyo  +  ^zo  -  (x-a)  F^  +  M^ 

a<x<  a  +b 

(B-10) 

X  =  a  +  b 

yo 


Byi  Fa’ 


(B-11) 


where  E  is  Young’s  modulus,  I  is  the  area  moment  of  inertia,  (M^q,  )  and  (Fy^,  Fyj )  are  the  moments 

and  forces  exerted  on  the  ends  of  the  beam. 


Three  possible  constraints  are  considered  at  each  end  of  the  beam.  At  the  front  end,  x  =  o,  these  are: 


and 


(a)  free  Fyo  =  ^  =  0 

(b)  pinned  yo  =  0 


(c)  fixed  y  =  0  ;  =  0. 

dx 


22 


Similar  conditions  can  be  assigned  to  the  rear  end,  x  =  a  +  b.  The  values  of  the  influence  coefficients  are 
given  in  Table  1  for  six  combinations  of  end  conditions;  fixed-fixed,  fixed-piimed,  pinned-fixed,  fixed- 
free,  free-fixed,  and  pinned-pinned.  The  construction  of  this  table  can  be  illustrated  by  considering  the 
rearward-facing  cantilever  case.  For  this  case,  Fyg  =  -F^,  =  -M^  -  aF^  (i.e.,  the  beam  is  fixed-free). 

Equation  (B-9)  can  be  integrated  twice  to  yield 


and 


El  .^  =  -  F 

dx  2 


+  xM^  o  <  X  <  a 


Fa  +  aM, 


a<  X  <a  +  b 


(B-12) 


At  -X^(X  -3a) 

Ely  = - 1- - i  Fj  + 


M, 


a^(3x -a) 
6 


P  +  a(2x  _a_) 

d. 


o<  X  <a 


a<  X  <  a  +  b. 


(B-13) 


At  X  =  a,  equations  (B-12)  and  (B-13)  become 

Ely(a)  =EI«ey  =  .^F,  + 


(B-14) 


and 


dx 


+  aM^. 


(B-15) 


The  above  coefficients  of  F^  and  are  precisely  those  given  by  the  fixed-free  entries  in  Table  1. 
The  other  coefficients  can  be  computed  in  a  similar  manner. 


23 


L 


Intentionally  left  blank. 
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LIST  OF  SYMBOLS 


FyC’  ^ZC 


(l/2)pS^Vpb[CM^  ^k^CN] 

(l/2)pS?V  [C„^  *  -  kfc^J 

(l/2)pSJv2c^_ 

distance  from  the  forward  end  of  the  beam  to  the  internal  component 
beam  influence  coefficients;  Table  1 

^xc  Pc  ~  ^tc  ^  1 
nic  Xc  « 

distance  from  the  rear  end  of  the  beam  to  the  internal  component 
beam  coefficients  derived  from  the  a^j’s. 

Magnus  moment  coefficient 

sum  of  the  damping  moment  coefficients 
static  moment  coefficient 
normal  force  coefficient 
center  of  mass 

nondimensional  beam  damping  orientation  coefficient,  equation  (19) 
nondimensional  beam  damping  displacement  coefficient,  equation  (18) 
(yc  +  iZc)/  * 

Young’s  modulus 

lateral  components  of  the  force  exerted  by  the  elastic  beam  on  the 
component 
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Ix’It 


^xb>  ^ 
^xc’  ^ 

K 


m 

“ifa 

P 
Pb 
Pc 
Q 
S 
s, 
t 

V' 


g 


axial  and  transverse  moments  of  inertia  of  the  total  projectile  relative  to  its  CM 
axial  and  transverse  moments  of  inertia  of  the  body  relative  to  its  CM 
axial  and  transverse  moments  of  inertia  of  the  component  relative  to  its  CM 

III 

/m^ ,  the  axial  radius  of  gyration 
,  the  transverse  raditis  of  gyration 


Pb  ^xb  Pc  ^xc 

reference  length 

lateral  components  of  the  moment  exerted  by  the  elastic  beam  on  the 
component 

mb  +  m^ 

mass  of  the  body 
mass  of  the  component 
projectile  spin  rate 
body  spin  rate 

component  spin  rate 

*  • 

0  +  i\|/ 

7c(P/4,  reference  area 

gyroscopic  stability  factor,  equation  (7) 

time 

projectile  velocity 

axial  distance,  measured  positive  rearward,  where  x  =  0  at  the  forward  end 
of  the  beam,  x  =  a  at  the  component,  and  x  =  a  +  b  at  the  rear  end  of  the 
beam 

distance  from  the  CM  of  the  (body  +  component)  to  the  CM  of  the  body 
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distance  from  the  CM  of  the  (body  +  component)  to  the  CM  of  the 
component,  where  m^x,,  +  m^.Xg  =  0 

transverse  displacements  of  the  component  CM  relative  to  the  body’s  axis 
of  symmetry 

the  projectile’s  angles  of  attack  and  sideslip  in  a  nonrolling  system 
Yy  +  i  Yz 

|r| 

orientation  angles  defining  the  component  axis  of  symmetry  with  respect  to  the 
projectile  axis  of  symmetry 

|E| 

projectile  pitch  rate 
♦ 

Kj/Kj,  a  damping  coefficient 

P  +  i  a 
air  density 

polar  angle  of  ^ 

fast  coning  fiequency 
r  phase  angle,  equation  (9) 

E  phase  angle,  equation  (10) 

projectile  yaw  rate 
natural  frequency  of  beam 


Intentionally  left  blank. 
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